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; ABSTRACT 

■ Context. Ion-neutral collisions may lead to the damping of Alfven waves in chromospheric and prominence plasmas. Neutral helium 
^— i " atoms enhance the damping in certain temperature interval, where the ratio of neutral helium and neutral hydrogen atoms is increased. 

| Therefore, the height-dependence of ionization degrees of hydrogen and helium may influence the damping rate of Alfven waves. 

Aims. We aim to study the effect of neutral helium in the damping of Alfven waves in stratified partially ionized plasma of the solar 
chromosphere. 

^■y^ Methods. We consider a magnetic flux tube, which is expanded up to 1000 km height and then becomes vertical due to merging with 

neighboring tubes, and study the dynamics of linear torsional Alfven waves in the presence of neutral hydrogen and neutral helium 
atoms. We start with three-fluid description of plasma and consequently derive single-fluid magnetohydrodynamic (MHD) equations 
<~j ■ for torsional Alfven waves. Thin flux tube approximation allows to obtain the dispersion relation of the waves in the lower part of 

[ tubes, while the spatial dependence of steady-state Alfven waves is governed by Bessel type equation in the upper part of tubes. 

Results. Consecutive derivation of single-fluid MHD equations results in a new Cowling diffusion coefficient in the presence of 
neutral helium which is different from previously used one. We found that shorter-period (< 5 s) torsional Alfven waves damp quickly 
in the chromospheric network due to ion-neutral collision. On the other hand, longer-period (> 5 s) waves do not reach the transition 
ffi ■ region as they become evanescent at lower heights in the network cores. 

Conclusions. Propagation of torsional Alfven waves through the chromosphere into the solar corona should be considered with 
caution: low-frequency waves are evanescent due to the stratification, while high-frequency waves are damped due to ion neutral 
collisions. 
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1 . Introduction al. I20071 De Pontieu et al. I20071 interested reader may find th e 

detailed review of oscillations in Zaqarashvili & Erdelyi 2009). 
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_i ■ Alfven waves play an important role in the dynamics of the Torsional Alfven waves do not lead to the displacement of 
\ solar atmosphere. Photospheric motions may excite the waves, magnetic tube axis, therefore they can be observed only with 
. . . which then propagate upwards along anchored magnetic field spectroscopic observ ations as periodic variation of spectral line 
and transport the energy into upper layers, where they may de- width (Zaqarashvili I2003l >. Observations of torsional Alfven 
posit the energy back leading to the chromospheric/coronal heat- waves wer e repo rted recently in c hromo spheric spectral lines 
ing and/or the acceleration of solar wind particles. Photospheric (Jess et al. I20091 De Pontieu et al. I2012I I. Upward propagating 
c3 \ magnetic field is concentrated in flux tubes, which are very dy- undamped Alfven waves can be also observed in the solar corona 
namic and may continuously change shape and/or merge due as an increase of non-thermal broad ening of coronal spectral 
to granular motions. Nevertheless, in a crude representation lines with height (Hassler et al. [1990l) . 
useful for theoretical modeling, the magnetic field can be ap- 
proximated as axis-symmetric tubes, therefore the excited pure The dynamics of Alfven waves in the chromosphere/corona 
Alfven waves are axis-symmetric i.e. they are torsional Alfven and their role in plasma heating are well studied (Hollweg |1981l 
waves. Considering a cylindrically symmetric flux tube, tor- [1984, Copil et al. 120081 Antolin and Shibata l2010l Vasheghani 
sional waves correspond to the azimuthal wavenumber set to Farahani et al. 120101 1201 II Morton et al. 12011V Thev can be ex- 
m — in the standard notation. On the other hand, granular cited due to the vortex motion at the photospheric level (Fedun 
buffeting may excite transverse magnetohydrodynamic (MHD) et al. 1201 laV In the case of inhomogeneous plasma Alfven 
kink waves in the tubes, which then can be transformed into lin- waves may also be excited due to resonances with other wave 
early polarized Alfven waves in upper layer due to the expansion modes (see, e.g., Soler et al. 120121 1. Torsional Alfven waves 
of tubes (Cranmer and van Ballegooiien 12005V The transverse can be used as supplementary tool for solar magnetoseismology 
oscillations have been observed with both, imaging and spec- (Zaqarashvili and Murawski 2007, Fedun et al. 1201 lbl Verth et 
troscopic observations (Kukhianidze et al. [2006, Zaqarashvili et al. l2010t . 
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In the solar chromosphere and prominences plasma is only 
partially ionized which leads to the damping of Alfven waves 
due to collision between ions and neutral hydrogen atoms (De 
Pontieu et al . 120011 Khodach enko et al. l2004l Leak e et al.|2005l 
Forteza et al [20071 Soler et al. |2009l Carbonell et al. l20T0l Singh 
and Krishan 120101) . Upward propagating Alfven waves may also 
drive spicules through ion-neutral collisions (Haerendel 119921 
James & Erdelvi l20U2"l Erdelyi & James [2004b . 

Single-fluid MHD description is a good approximation for 
low-frequency waves, but it fails when the wave frequency ap- 
proaches to ion-neutral collision frequency, and consequently 
the multi-fluid MHD description should be used. It was shown 
by Zaqarashvili et al. (1201 1 al l that the damping rate is maximal 
for the waves those frequency is near ion-neutral collision fre- 
quency i.e. higher and lower frequency waves have less damping 
rates. 

Beside the neutral hydrogen atoms, the solar plasma may 
contain significant amount of neutral helium atoms, which may 
enhance the damping of Alfven waves. Soler et al. ( 120101 ) sug- 
gested that the neutral helium has no significant influence on 
the damping rate in the prominence cores with 8000 K temper- 
ature. On the other hand, Zaqarashvili et al. ( 1201 lbl l showed 
that the neutral helium may significantly enhance the damping 
of Alfven waves in the temperature interval of 10000-40000 K, 
where the ratio of neutral helium and neutral hydrogen atoms is 
increased. This means that the helium atoms can be important in 
upper chromosphere, spicules and prominence-corona transition 
regions. The calculation of neutral helium effects in Zaqarashvili 
et al. (1201 lbl l was done in homogeneous medium with uniform 
magnetic field and ionization degree. On the other hand, the 
ion and neutral atom number densities and consequently ioniza- 
tion degree are significantly changed with height in the chromo- 
sphere, which may influence the damping of Alfven waves due 
to ion-neutral collisions. 

Stratification may significantly influence the dynamics of 
waves leading to their reflection at transition region, wave mu- 
tual transformation and/or evanescence. It introduces a cut-off 
frequency for Alfven waves in fully ionized plasma, i.e. waves 
with lower frequency than the cut-off value are evanescent in the 
solar atmosphere (Musielak et al. 119951 Murawski & Musielak 
l20T0l >. 

In this paper, we study the effects of stratification and neu- 
tral helium atoms on the propagation/damping of Alfven waves 
in the solar chromosphere. We consider a magnetic flux tube, 
which is expanded up to 1000 km height and then becomes ver- 
tical due to merging with neighboring tubes. Consequently, we 
consider torsional Alfven waves, which are only pure incom- 
pressible waves in tubes. 



2. Linear torsional Alfven waves 

We consider a vertical magnetic flux tube embedded in the strat- 
ified solar atmosphere. We use cylindrical coordinate system 
(r, 6, z) and suppose that the unperturbed magnetic field, Bq, is 
untwisted i.e. Bog = 0. Plasma consists in electrons, protons 
(H + ), singly ionized helium (He + ), neutral hydrogen (H) and 
neutral helium (He) atoms. We consider the Alfven waves po- 
larized in 9 direction i.e. the only non-zero components of the 
perturbations are vg and be- Therefore, the waves are torsional 
Alfven waves. We start with the three-fluid approach for par- 
tially ionized hydrogen-helium plasma, when one component is 
ion+electron gas and other two components are neutral hydro- 
gen and neutral helium gasses (Zaqarashvili et al. 1201 lbl l. Then, 
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Fig. 1. Height dependence of atmospheric parameters according 
to FAL93-F model (Fontenla et al. 1993): magenta, red and blue 
lines correspond to proton, neutral helium and neutral hydro- 
gen number densities, respectively. Green line is the total (pro- 
ton+neutral hydrogen+neutral helium) number density. 

the linear Alfven waves are governed by the equations (magnetic 
diffusion due to ion-electron collision is neglected): 

d, B a -V a H + a He a R a He , W1 , 

-^-(rve) = (rbg) (rv e )+ — (rv m )+ (rv He e),(l) 

dt 4-npi pi pi pi 

-r(rv m ) = (rvg) (rv m ) + ("'Hee), (2) 

ot pho Pho Pho 

d . . a H , a He + a He H , , «HeH . . 
-^(rvHefl) = ^ rVe > Ov H ee) + (rv m ), (3) 

Ot PHeO PHeO PHeO 



£-«*-v,(2). 



(4) 



where Vg (vne, vn e e) are the perturbations of ion (neutral hydro- 
gen, neutral helium) velocity, bg is the perturbation of the mag- 
netic field, and p, (pm, PHeo) is the unperturbed ion (neutral hy- 
drogen, neutral helium) density. Here we use the definitions 



CH = <*H+H + CHe+H, ^He = CH+He + O'He+He, 



(5) 



where a denotes the coefficient of friction between different sort 
of species. The collision between electrons and neutral atoms is 
also neglected as it has much less effect in the damping of MHD 
waves than the collision between ions and neutrals. In the present 
model the motions of torsional Alfven waves are normal to the 
direction of gravity, so that gravity does not explicitly appear in 
the equations. The effect of gravity is indirectly present in the 
equations through the stratification of plasma density. 

The coefficient of friction between ions and neutrals (in the 
case of same temperature) is calculated as (Braginskii ll965l ) 



nin n m in a in -, 



8kT 



(6) 



where T is the plasma temperature, m, (m„) the ion (neutral 
atom) mass, m,„ = m/m„/(m,- + m n ) is reduced mass, n,- (n„) is 
the ion (neutral atom) number density, cr,„ = 7r(r,- + r„) 2 » nr\ 
is ion-neutral collision cross section for elastic collisions from 
Braginskii d!965l l (mean atomic cross section is nr\ = 8.7974 x 
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1(T 17 cm 2 ), and k = 1.38 x 1(T 16 erg K -1 is the Boltzmann con- 
stant. 

From here we adopt the single-fluid approximation. To jus- 
tify whether or not the single-fluid approximation is valid in the 
solar atmosphere let us estimate the ion-neutral collision fre- 
quency. Mean collision frequency between ions and neutrals can 
be calculated as (Zaqarashvili et al. 1201 lbb 



Vin = Ctb 



1 1 
1 



mini m„n, 



(7) 



The collision frequency is very high in the photosphere and de- 
creases upwards. For example, The collision frequency between 
protons and neutral hydrogen can be estimated as 8.6 • 10 6 Hz, 
6.2 • 10 3 Hz and 24 Hz, at z = 0, z = 900 km and z = 1900 km 
correspondingly (ion and neutral atom number densities were 
taken from Fontenla et al. 119931 model-F). Therefore, the Alfven 
waves with periods > 1 s can be studied in the single-fluid ap- 
proach. 

To obtain the governing equations in the single fluid approx- 
imation, we consider the total density 



P - Pi+PH + PHe 

and velocity of center of mass 

PiVg + PHVm + PHeVRee 



Vg = 



Pi +PH + PHe 



(8) 



(9) 



We also consider the relative velocity of ion and neutral hydro- 
gen as whs = vg - Vhs and the relative velocity of ion and neutral 
helium as wn e g = vg - VH e e- Then one may find that 



Vg = Vg + % H Wm + fee-VVn 



(10) 



where % H = p H /p and g He = p He /p- 

Consecutive substraction of Eq. (Q]) and Eq. (0, Eq. (0 and 
Eq. (0, Eq. (0 and Eq. (0 and neglecting the inertial terms for 
relative velocities lead to the equations 



rw m = 



rWHe 



a H e a HeH (i . . 
Sfl + (?fl + SHe) 

a a 



a H , <XHeH lt . - . 
SHe + (?H + SHe) 

a a 



go-V 
An 

4-7T 



(rbg), 



<rb e ), 



(11) 



(12) 



where a = anotue + anaHeH + OHe^HeH- The neglecting of the 
inertial terms is justified for waves those frequency is less than 
ion-neutral collision frequency. This is the key step to transform 
the multi-fluid equations into the single fluid description. 

Then the sum of Eqs. ([TJ-OJ and using Eqs. ©, ( fTOt . (ITTb 
and ([T2I 1 lead to the single-fluid equations 



^ =r(Bo .V)(^) + r(B -V) 



ric(z) B V 



(rbg) 



(14) 



where 



4k a H a He + a H a He H + a He UHeH 



is the coefficient of Cowling diffusion. The first two terms in 
the nominator of Eq. ( TT3T > are due the frictions of neutral hydro- 
gen and neutral helium separately with regards to ions, while the 
last term is the friction of neutral hydrogen+helium fluid with 
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Fig. 2. The ratio of Cowling diffusion coefficient and Alfven 
speed square, iJcjVl, vs height according to FAL93-F model 
(Fontenla et al. [T993l . 

regards to ions. Therefore, the last term is important in weakly 
ionized plasmas, while the first two terms are important for rel- 
atively high value of ionization degree. Soler et al. (120101 used 
the following expression for Cowling diffusion in the presence 
of neutral helium 

B 2 Q {tH+tHef .... 
Tic = -. 7 (16) 

47T an + a He 

which can be obtained from our expression (Eq. \T5[ when 
oh, aHe ^ &HeH i.e. when plasma is only weakly ionized. This is 
so because Soler et al. ((2010) considered in their derivation that 
both neutral hydrogen and neutral helium have the same velocity. 
That is equivalent to take au e n 00 ■ Therefore the present de- 
scription is more general than that of Soler et al. (1201 Oi l. The ex- 
pression used by Soler et al. ( 1201 Oi l is only approximately valid 
for the photosphere and lower chromosphere, while in the upper 
chromosphere, spicules and prominences the general expression 
Eq. (fT~5T > should be used. 

Let us now consider a new coordinate s, which represents 
a longitudinal coordinate along the unperturbed magnetic field, 
Bo, then Eqs. (fT3l- (fT4l i can be combined into the single equation 
in terms of this new coordinate, namely 



B 0s d 



(13) Anpr 2 ds 



8 2 Ug 
dt 2 



r Bot- 
es 



Bo 



Anpr 2 ds 
'4npT] c dU g 



r 2 B Qs 



dUg 



Bi 



dt 



(17) 



where Ug = Vg/r. Note that r is a function of s. This equa- 
tion is simplified near the axis of symmetry, where Bo s (s)r 2 (s) w 
constant (Hollweg |19811 l, which is the conservation of magnetic 
flux near the tube axis, where field lines are almost perpendicu- 
lar to the tube cross-section. Then this equation may be written 
as 



(15) d 2 Ug 



3t 2 
where 

V A = 



1 + 



Jlc(s) d 
V 2 A (s) dt 



Bqs(s) 
y/4np(s) 



(18) 



(19) 
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Fig. 3. Vertical magnetic flux tube model: below 1000 km the 
Alfven speed is constant due to the thin flux tube approximation 
and above 1000 km height the Alfven speed increases exponen- 
tially. 



is the Alfven speed. 

Solar atmospheric parameters depend on the height due to 
gravity. Plasma ionization degree is also a function of altitude 
due the increase of temperature with height. The plasma is only 
weakly ionized in the photosphere/lower chromosphere. The in- 
crease in temperature with height leads to the ionization of hy- 
drogen and helium atoms, which become almost fully ionized in 
the solar corona, so the transition occurs near the region of sharp 
temperature rise. 

Figure 1 shows the density of different species vs height 
according to the FAL93-F model (Fontenla et al. 119931 ). This 
model includes the dependence of ionization degree on heights 
for both hydrogen and helium. The neutral atom number densi- 
ties are much higher than the ion number densities at the lower 
heights, but become comparable near ~ 1900 km, which corre- 
sponds to the temperature of about 10000 K. It is seen that the 
neutral atoms are more stratified than protons between 500 and 
2000 km heights: the neutral hydrogen, neutral helium and to- 
tal number densities have scale height about 180 km, while the 
scale height of protons is much larger. The total number density 
is similar to the neutral hydrogen number density until ~1500 
km. Above this height they start to diverge due to the increase of 
ionization degree. 

Fourier analysis of Eq. ( TT8l with exp[-/a;?], where u> is the 
wave frequency, gives the equation 



ds 2 



1 vfr) 



L0 2 Uff 



(20) 



which governs the spatial dependence of steady-state torsional 
Alfven waves. The type of the equation depends on the depen- 
dence of Va and r\ c on s. 

Very important parameter in Eq. d20T i is the ratio n c /V^, 
which does not depend on the magnetic field structure with 
height as B 2 )s is canceled. We plot the ratio vs height according 
to VAL93-F model on Fig. 2. It is clearly seen that the ratio is al- 
most constant along whole chromosphere and abruptly increases 
near the transition region. Therefore, in order to find analytical 



solutions to Eq. (Got we consider the ratio as a constant along s 
in the chromosphere 



const. 



(21) 



From Figure 2 we see that this approximation seems to be 
good in the upper part of the chromosphere, i.e., between 800- 
2000 km, while it may be not so accurate in the lower part 
of the chromosphere, i.e., in the interval between 500-800 km. 
Nevertheless, we use this approximation along whole chromo- 
sphere for convenience. 

Plasma B (-d,np/B^) is constant with height in the isother- 
mal atmosphere in thin vertical magnetic flux tubes, when tem- 
peratures inside and outside the tubes are same (Roberts |2QQ4b . 
In this case, the Alfven speed is also constant with height yield- 
ing Va(s) = Vao = const. The constancy of the Alfven speed 
in thin tubes is caused due to the compensation of density vari- 
ation by the magnetic field strength variation with height. The 
thin flux tube approximation is valid up to 1000 km above the 
photosphere (Hasan et al. 120031 Cranmer and van Ballegooijen 
120051 1. Above this height, the magnetic tubes are thick and they 
probably merge with neighboring tubes producing almost ver- 
tical magnetic field (Fig. 3). Then, the Alfven speed may vary 
with height due to the decrease of density, while the magnetic 
field is constant. According to this model, we consider the verti- 
cal magnetic tube which consists in two different parts: the lower 
part with constant Alfven speed and the upper part with expo- 
nentially increasing Alfven speed. Then we study the torsional 
Alfven waves in the two parts of magnetic flux tube separately. 



3. Lower chromosphere: thin flux tube 
approximation 

In the lower part of magnetic flux tube, the thin tube approxima- 
tion is used which yields Va(s) = const (see previous section). 
Note that the thin tube approximation is still compatible with a 
magnetic flux tube expanding with height. The only restriction 
is that the wavelength remains much longer than the tube radius. 
This case was also studied by Soler et al. J2009) for the case of 
torsional Alfven waves in prominence threads. The constancy of 
Alfven speed means that the coefficient of Cowling diffusion is 
also constant. Then the coefficients of Eq. (1201 do not depend 
on s, therefore Fourier analysis with exp[i(k s s - cot)] gives the 
dispersion relation 



co 2 + irj c k 2 u> 



■ k 2 V 2 
K i v A 



0. 



which has two complex solutions 



+k s V A 



tl 2 kg _ . Vck s 



(22) 



(23) 



The dispersion relation and its solutions agree with Equations 
(33) and (34) of Soler et al. (120091) . Real part of this expression 
gives the wave frequency, which shows that there is the cut-off 
wavenumber due to Cowling's diffusion, k s = 2Va/t] c . The cut- 
off wave number due to ion-neutral collision was found in differ- 
ent situations of the solar prominences (Forteza et al. 120071 Soler 
et al. 120091 Barcelo et al. 1201 U . It is the result of approximation 
when one transforms multi-fluid equations into the single-fluid 
MHD and therefore it has no physical ground (Zaqarashvili et 
al. l20T2l l. 
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On the other hand, the imaginary part of this expression gives 
the normalized damping rate as 



m\ = 



OJ, 



k s V A 



1 k s V A a H eP 2 H + a H p 2 He + CtHeHiPH + PHef 



P 



a H a He + a H a He H + a He a Heli 



,(24) 



where we have used the full expression of Cowling's coefficient, 
rj c . The plasma is almost neutral in the low chromosphere, there- 
fore (%HeH ^ atH, ctHe- This leads to the normalized damping rate 



1 k s V A 

2 p 



(PH + pHe) 



a H + a He 



(25) 



This expression was used by Soler et al. (2010) to calculate the 
effect of neutral helium in prominence cores. 

On the other hand, when the number density of neutral atoms 
is much less than ion number density, then ctHeH <K <*h, <iHe and 
the normalized damping rate is 



1 k s V A 



2 2 

an a He 



(26) 



When ion number density is of the same order as the num- 
ber density of neutral atoms (like in spicules or in prominence 
cores), then the general expression of damping rate (Eq. \2M 
should be used. The expression Eq. [24]may significantly change 
the contribution of neutral helium atoms. 

Fig. 4 shows the normalized damping time {T^/To = l/d>,-) 
vs the period of Alfven waves (Tq = 2n/k s V A ) for three 
different areas: a) FAL93-A corresponding to faint cell cen- 
ter area, b) FAL93-F corresponding to bright area of chromo- 
spheric network, and c) prominence cores. We use the follow- 
ing parameters: a) n,=3.26xl0 10 cm 4 , «h=3.01x10 13 cm -3 , 
n He =3.01xl0 12 cm 3 taken from FAL93-A at 975 km height 
corresponding to as 5480 K temperature; b) «,-=2.49xlO n cm -3 , 
«h=5.24x10 13 cm -3 , «H e =5.26xl0 12 cm -3 was taken from 
FAL93-F at 975 km height corresponding to 6180 K temperature 



and c) n,=10 cm , « w =2xl0 cm , « We =2xlO cm" 3 was 
taken for prominence cores corresponding to 8000 K tempera- 
ture. It is clearly seen that the presence of neutral helium atoms 
significantly enhances the damping of Alfven waves. It is also 
seen that the damping of torsional Alfven waves is more efficient 
in bright cell center area in the chromosphere and in prominence 
cores, where the waves are damped over few wave periods. On 
the other hand, the damping of torsional Alfven waves due to 
ion-neutral collision is inefficient in bright chromospheric net- 
work center. The clear difference of damping in these two areas 
can be understood in terms of ion-neutral collision frequency. In 
the absence of neutral helium atoms, Eq.|24]can be rewritten as 
(Zaqarashvili et al. l2011bb 



1 k s V A pH 

\0>i\ « , 

2 V iH pi 



(27) 



which clearly indicates that the normalized damping rate de- 
pends on the ratio of Alfven frequency over ion-neutral hydro- 
gen collision frequency (v,h) and the ratio of neutral and ion 
fluid densities. For p, «; p#, which is the case at 975 km height, 
ion-neutral hydrogen collision frequency is proportional to neu- 
tral hydrogen density, vm ~ pu- This means that the normalized 
damping rate is inversely proportional to ion number density, 
Pi. Proton number density is almost 10 times smaller in the cell 
center than in the network center, which leads to almost 10 times 
difference in damping (Fig. 4). 

In the first part of this section we have studied the temporal 
damping of Alfven waves, i.e., we have set a real k s and solved 



the dispersion relation to obtain the complex a>. In the follow- 
ing paragraphs we study the spatial damping. This means that 
the dispersion relation (Eq. l22l is now solved for the complex 
wavenumber, k s , assuming a fixed and real a>. The solution of 
the dispersion relation for the square of k s is 



kt = 



of 



V 2 - iUT] c 



(28) 



We use this expression to find the real and imaginary parts of k s , 
assuming k s — k r + ik t . After some algebraic manipulations we 
can obtain the expressions for kf and ki, namely 



2t/2 



oj z V- 



2t/2 



ofV 



2(y4 + oj 2 tf) 2(V A + coWc) 



1 + 



V A 



ki = 



1 



VcOJ 



k r 2(V.+coWc) 



(29) 



(30) 



The expression of k r is computed from Eq. ( 129b as k r — ± -y/fc 2 . 
Note that the + sign in front of the second term in the expression 
of k 2 r leads to 4 possible values of k r . The two values of k r with 
+ sign in the expression of kf. correspond to propagating waves, 
namely 



k r = ± 



\2{V\+oj^I) ' 2(y4 + w 2^ } 



1 + 



(31) 



where the + sign in front of the square root refers to upward and 
downward waves, respectively. 

To obtain more simplified analytical formulas let us consider 

the case of weak damping, i.e., ^f- <K 1. The expression of k r 

A 

for propagating waves becomes 



k r 



Vl + <oW 



1 + 



8 V? , 



(32) 



We now take the + sign corresponding to upward propagating 
waves and use this expression to compute the value of A:, in the 

limit <k 1, namely 

'A 



ki 



1 Tlcor 
2k r V A 



2V\ 



(33) 



The corresponding normalized damping rate using the full ex- 
pression of rj c is 



1 oj a He p 2 H + a H p 2 He + a HeH (p H + p He ) 2 



(34) 

OJ/ V A 2 p a H a He + a H a He H + a He a HeH 

which is equivalent to Equation (124-b obtained for the temporal 
damping if we replace co by k s V A . Therefore, both spatial and 
temporal damping rates of the Alfven waves are equivalent. 

On the other hand, the two values of k r corresponding to 
the - sign in the expression of kf. (Equation (l29l ) correspond 
to evanescent waves, namely 



k r 



+ 0jWc) 2(^+0,^2) 



1 + 



K > 



(35) 



which for weak damping, i.e., 2iLj- <g l ; becomes 



k T = ±i 



orrf c 



(36) 
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This corresponds to evanescent perturbations that do not prop- 
agate from the location of the excitation, and so these solutions 
are not relevant for the present study. 

4. Upper chromosphere: exponentially increasing 
Alfven speed 

In the upper part of magnetic flux tube, where the neighboring 
tubes are merged, we assume that the Alfven speed is exponen- 
tially increasing (note that the ratio of Cowling diffusion coeffi- 
cient and Alfven speed square is assumed to be constant again) 



V A (s) = V AO exp^j, 



(37) 



where h is the scale height and V^o corresponds to the value of 
Alfven speed at the height of 1000 km (i.e. s — also corre- 
sponds to the height of 1000 km). In this case, Eq. ( 1201 can be 
rewritten as 



d 2 U e 

ds 2 

where 



k 2 

K s0 



+ exp 



o. 



or 



V 2 



(38) 



(39) 



A0 - 10JT] M 



and rjcQ is the value of Cowling diffusion at the height of 1000 
km. The solution of this equation is 



Ug = Ug(0)F (O, 



(40) 



where Fq is Bessel, modified Bessel or Hankel functions of zero 
order and 



f — 2hk S Q exp 



s 

'2h 



Propagating waves are described by Hankel functions 
(H^iO = J (O + iY (O, H 2 Q (0 = MO - iY (O) and stand- 
ing waves are described by Bessel functions (Jq(^)). Hn(€) an d 
Hq(() functions have the same spatial dependence in the case of 
fully ionized plasma (Fig. 5). The upward and downward propa- 
gating waves can be distinguished by the sign of time-averaged 
Pointing flux (Hollweg 11984b . which shows that H 2 ^) is up- 
ward propagating and Hq(0 is downward propagating waves. In 
the case of partially ionized plasma distinguishing between up- 
ward and downward propagating waves is complicated due to 
the complex argument of Hankel functions (De Pontieu et al. 
1200 U . However, it is still possible to distinguish them by care- 
ful look into the height-structure of waves (Fig. 5). Red line, 
which corresponds to H^((), has smaller amplitude at higher 
heights comparing to the fully ionized case (green line), while 
blue line, which corresponds to Hq(£), has the stronger ampli- 
tude. This means that //q(^) governs the wave, which damps at 
higher heights due to ion-neutral collisions, therefore it governs 
the dynamics of upward propagating waves. On the same basis, 
Hq(() governs the downward propagating waves. The standing 
waves, expressed by Jq{0, can be formed after superposition of 
upward and downward propagating waves. 

A mathematical proof that the function corresponds to 
upward waves is obtained from the asymptotic expansion of the 
Hankel functions for large arguments (see, e.g., Stenuit et al. 
1 1999b . For large £, which eventually means small s, the depen- 
dence of Hankel functions on s can be written as 

HliO ~ expHfcok), Hl(0 ~ exp(i\k sQ \s), 



where k s o is given in Equation d39b . These expansions clearly 
show that corresponds to upward propagating waves when 
the temporal dependence is exp[-;wf]. In the rest of the paper 
we consider only upward propagating waves. Note that only the 
real parts of Hl(Q and F&(£), i.e., the part of the solution with 
the physical meaning, are shown in all plots. 

Fig. 6 shows the height-dependence of upward propagat- 
ing torsional Alfven waves with different periods in the bright 
chromospheric network cores. At 1000 km height we use the 
number densities from FAL93-F model as n,=3.22xlO n cm -3 , 
«#=2.65xl0 13 cm -3 , «H e =2.68xl0 12 cm 4 . The plasma temper- 
ature is set as r=8000 K and the Alfven speed is set as Vao= 
20 km s . It is seen that the behavior of longer period waves 
( > 10 s) is not significantly affected by ion-neutral collisions: 
the height dependence is similar for fully ionized (green line) 
and partially ionized (blue and red lines) plasmas (upper panel). 
The shorter period waves display more significant dependence 
on ion-neutral collision frequency. For example, the torsional 
Alfven waves with period of 3 s are significantly damped in par- 
tially ionized plasma comparing to the fully ionized case (lower 
panel). It is also seen that the presence of neutral helium (red 
line) enhances the damping of Alfven waves compared to neutral 
hydrogen only (blue line). Fig. 6 shows that the spatial depen- 
dence of waves are not oscillatory above particular heights i.e. 
waves become evanescent. In order to study this phenomenon, 
we consider wave propagation in fully ionized plasma. 

Fig. 7 displays the velocity of torsional Alfven waves vs 
height for different wave periods. It is seen that torsional Alfven 
waves with period of 20 s become evanescent above 1500 km 
height. Waves with 5 s period become evanescent above 1800 
km height, but the waves with 1 s period penetrate up to 2600 km 
height. Therefore, the longer period waves ( > 5 s) do not reach 
to the transition region (2000 km height), but become evanes- 
cent at lower heights. This gravitational cut-off of Alfven waves 
is known for a long time (Musielak and Moore, 1995 ). Recently, 
Murawski and Musielak (2010) also showed that the linearly po- 
larized Alfven waves become evanescent above some heights in 
fully ionized plasma. Our result confirms their analyses. 

5. Discussion 

Propagation of Alfven waves in the chromosphere is very im- 
portant for chromospheric/coronal heating as they may carry the 
photospheric energy into the upper layers. The magnetic field 
is concentrated in flux tubes at the photospheric level, there- 
fore the torsional Alfven waves can be generated due to vortex 
motions (Fedun et al. 1201 lab . Observation of torsional Alfven 
waves is possible through the periodic variation of spectral line 
width (Zaqarashvili 2003). Recently, Jess et al. (2009) reported 
the observation of torsional Alfven waves in the lower chromo- 
sphere. The energy flux carried by the waves was enough to heat 
the solar corona. Therefore, it is interesting to study whether the 
waves may penetrate into the corona. 

Here we study the torsional Alfven waves in chromospheric 
magnetic flux tubes with partially ionized plasma taking into ac- 
count neutral hydrogen and neutral helium atoms. We consider 
the stratification due to gravity, therefore the tubes are expanding 
with height. Hasan et al. d20031 > showed that the magnetic flux 
tubes can be considered as thin up to 1000 km height from the 
surface. Above this height, the neighboring magnetic flux tubes 
merge and the field lines are almost parallel. We use this model 
and split the magnetic tubes into two layers. We use the expand- 
ing thin flux tube approximation in the lower chromosphere be- 
low 1000 km. In the upper layer, above 1000 km we consider the 
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magnetic field lines as vertical, so only medium density is chang- 
ing with height. We use FAL93 model (Fontenla et al. 11993b 
for the dependence of ion, neutral hydrogen and neutral helium 
number densities with height. 

We start with three-fluid MHD description of partially ion- 
ized plasma, where one component is electron-proton-singly 
ionized helium and the two other components are neutral hydro- 
gen and neutral helium gasses. Then we proceed to the single- 
fluid description by standard procedures and neglect the inertial 
terms of ion-neutral hydrogen and ion-neutral helium relative 
velocities. We obtain the Cowling diffusion coefficient (Eq.[T5ll. 
which is different than as previously used one. Namely, the ex- 
pression Eq. ( TT5b contains additional terms (the first and second 
terms in nominator) expressing the collision of ions with neu- 
tral hydrogen and neutral helium atoms separately. The previ- 
ously used expression of Cowling diffusion coefficient is valid 
in weakly ionized plasma, where the coupling between neutral 
atoms are stronger. 

In the lower part of magnetic flux tube, where the thin tube 
approximations is used, the dispersion relation for the torsional 
Alfven waves is obtained. The solution of the dispersion relation 
shows that the damping of Alfven waves due to ion-neutral col- 
lision is significant in faint cell center areas at all frequencies, 
while only high-frequency waves are damped in the chromo- 
spheric network cores (Fig. 4). This means that torsional Alfven 
waves may propagate without problem up to 1000 km above the 
surface in chromospheric network cores. The effect of neutral 
helium atoms is important in all cases and it increases for lower 
frequency waves. 

In the upper part of magnetic flux tube, the height depen- 
dence of torsional Alfven wave velocity is governed by Bessel- 
type equation, therefore the solutions are Bessel or Hankel func- 
tions with zero order and complex arguments. We consider only 
upward propagating waves, which are expressed by Hankel func- 
tions, Hq. The solution of steady-state torsional Alfven waves 
shows that the long-period waves are not affected by ion-neutral 
collisions. Only short-period waves (< 5 s) are damped signif- 
icantly (Fig. 6). Existence of neutral helium atoms enhances 
damping of the waves but not significantly. On the other hand, 
long-period waves show no oscillatory behavior above cer- 
tain heights, so they become evanescent in upper part of the 
chromosphere. This phenomenon is called gravitational cut-off 
(Musielak et al. 119951 Murawski & Musielak 120101 ). It is seen 
that only the waves with very short period (~ 1 s) may reach 
the transition region in the case of fully ionized plasma (Fig. 7). 
This means that the low-frequency waves can not reach the tran- 
sition region due to the gravitational cut-off, but high-frequency 
waves are damped due to ion-neutral collision. Then the tor- 
sional Alfven waves can not reach the transition region at all. 
This is true near the axis of magnetic flux tubes, where our ap- 
proximations are valid. However, it is possible that some wave 
tunneling exists there and the evanescent tail of the waves can 
indeed reach the transition region and the corona. Due to the 
change of the ambient conditions the waves may become propa- 
gating again in the corona. In the outer part of tubes, where mag- 
netic field lines are more inclined, the problem of wave propa- 
gation along the chromosphere may not arise at all. It requires 
further detailed study. 

On the other hand, Alfven waves may easily penetrate into 
the corona if they are excited in the higher part of the chro- 
mosphere (Hansen and Cally I2012I I. Long-period acoustic os- 
cillations (p-modes) are evanescent in the lower atmosphere due 
to the gravitational stratification, however they may propagate 
with an angle about the vertical and may trigger Alfven waves in 



the upper chromosphere through mode conversion (Zaqarashvili 
and Roberts 120061 Cally and Goossens 120081 Cally and Hansen 
1201 11 Khomenko and Cally 120121 . If low-frequency torsional 
Alfven waves are excited above 1000 km height, then our model 
also allows the propagation of the waves into the corona as the 
evanescent point can be shifted up above the transition region 
(Fig. 7). 

Recently, Vranjes et al. (J20081 suggested that the energy flux 
of Alfven waves excited in the photosphere is overestimated and 
thus the Alfven waves can be hardly excited in the photosphere 
(but see Tsap et al. 1201 ll for alternative point of view). This is in- 
deed an interesting question and should be addressed adequately. 
One possible explanation is that Vranjes et al. (120081) considered 
the magnetic field strength of 100 G in the photosphere, which is 
too weak for flux tubes, where the magnetic field strength is of 
order kG. Therefore, it is possible the Alfven waves are excited 
inside flux tubes only, but not outside the tubes, where the mag- 
netic field is relatively weak. This point needs further discussion. 

Ion-neutral collisions are important only for short-period 
waves (< 5 s). On the other hand, only long period transverse 
oscillations (> 3 min) are frequently observed so far in the so- 
lar atmosphere. However, some observations from ground based 
coronagraphs show the oscillations of spicule axis with periods 
of < 1 min (Zaqarashvili et al. [2007, Zaqarashvili & Erdelyi 
120091 ). Several reasons could be responsible for the absence of 
high-frequency oscillations: not enough tempo-spatial resolution 
of observations, hard excitation in the photospheric level, quick 
damping in the lower atmosphere etc. We hope that more so- 
phisticated observations will help us to improve our knowledge 
about high-frequency oscillations in the future. 

6. Conclusions 

Propagation of torsional Alfven waves along expanding vertical 
magnetic flux tubes in the partially ionized solar chromosphere 
is studied taking into account the ion collisions with neutral hy- 
drogen and neutral helium atoms. The waves propagate freely 
in the lower chromosphere up to 1000 km and may damp by 
ion-neutral collisions. The new expression of the Cowling dif- 
fusion (and consequently damping rate) including the neutral 
helium is obtained (Eq. 15). Neutral helium atoms may sig- 
nificantly enhance the damping of Alfven waves. On the other 
hand, the long period (> 5 s) waves propagating near the tube 
axis become evanescent above some height in the upper chro- 
mosphere not reaching the transition region, while the short pe- 
riod waves (< 5 s) are damped by ion-neutral collisions. Hence 
the torsional Alfven waves may not reach the transition region 
and consequently the solar corona unless some tunneling effects 
are considered. This means that the coronal heating by photo- 
spheric Alfven waves should be considered by caution. At the 
same time, energy of the waves can be dissipated in the chromo- 
sphere leading to the heating of ambient plasma. 
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Fig. 4. Normalized damping time of torsional Alfven waves 
(Td/To) vs wave period (To) for three different situations from 
top to bottom: faint cell center area (FAL93-A), bright network 
(FAL93-F) and prominence core. The blue line indicates the 
wave damping due to collision of ions with neutral hydrogen 
atoms only. While the red line indicates the wave damping due 
to collision of ions with neutral hydrogen and neutral helium 
atoms. 
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Fig. 5. Height dependence of steady state torsional Alfven waves 
with period of 5 s. Green line corresponds to Hq(£) and Hq(() 
from Eq. ( f40b in fully ionized plasma. Blue and Red lines cor- 
respond to Hq(() and Hq(£), respectively, in partially ionized 
plasma. 
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Fig. 6. Height dependence of steady state upward propagating 
torsional Alfven waves with periods of 20 s (upper panel), 10 
s (middle panel) and 3 s (lower panel). Green lines correspond 
to the fully ionized plasma, blue (red) lines correspond to par- 
tially ionized plasma with neutral hydrogen atoms (neutral hy- 
drogen+neutral helium atoms). 
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Fig. 7. Height dependence of steady state upward propagating 
torsional Alfven waves with periods of 20 s (red line), 5 s (blue 
line) and 1 s (green line) of different periods in fully ionized 
plasma. 



